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A quantum kinetic theory for correlated charged— particle systems in strong time- 
dependent electromagnetic fields is developed. Our approach is based on a sys- 
tematic gauge— invariant nonequilibrium Green's functions formulation. We con- 
centrate on the selfconsistent treatment of dynamical screening and electromag- 
netic fields which is applicable to arbitrary nonequilibrium situations. Numerical 
results for the nonlinear plasma heating by the laser field and the electron— ion col- 
lision frequency including multi-photon absorption (inverse bremsstrahlung) are 
presented. 



C , 1 Introduction 

With the progress in short -pulse laser technology, high intensity electromag- 
■ netic fields are now allowing to create strongly correlated quantum plasmas 

Q\ \ in extreme nonequilibrium conditions with unprecedented applications . At 

the same time, optical techniques for time-resolved diagnostics are improving 
remarkably || creating the need for a quantum kinetic theory of dense non- 
ideal plasmas in intense laser fields. A central problem is the simultaneous 
account of correlations (scattering) between the particles and of the influence 
of the electromagnetic field. On the other hand, such a theory is of interest 
Q for quantum transport phenomena in solids subject to intense THz fields, e.g. 

where the common rotating wave approximation cannot be applied. 

In this paper we take advantage of the Kadanoff-Baym formalism to treat 
charged particles and plasmons in a fully symmetric way. The time-dependent 
electromagnetic field is treated classically. Using a gauge-invariant extension 
of the generalized Kadanoff-Baym ansatz, we obtain a closed equation for the 
Wigner distribution which is solved for various limiting cases. 
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2 Quantum kinetic equations for charged particles and longitudinal 
photons 



We consider charged particles (electrons and ions/holes/positrons etc.) inter- 
acting via the longitudinal Coulomb force which is equivalently described in 
terms of emission and absorption of longitudinal photons. The nonequilibrium 
state of this system is conveniently characterized by correlation functions of 
fermions of species a, go (1, 1') and longitudinal photons (screened potential) 
M,V S <(1,1'), which are defined as averages of the respective field operator 
pairs [^), where we use the notation 1 = Ti,tx. In particular, the density 
matrices of carriers and longitudinal photons follow from the " < " functions, 



/a(r 1 ,r' lJ *i) = -tfiff<(l J l / )|t 1 =t i andiV(r 1) r' 1 ,t 1 ) = -ifi7 s< (l ) l')| il = ti . The 



transverse electromagnetic field is given by the vector potential A and will be 
treated classically, i.e. it obeys Maxwell's equations (||). For completeness, we 
include also the longitudinal field which is due to external sources, cxt , Q. 
The time evolution of the correlation functions is determined by the Kada- 



noff-Baym equations Rj 




dt\ 2m a \ i 
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A^(M')=£ 



dl nf c (i,T)T/|(I,i')-n| c (i,i)Kf (i,i') ,(2) 
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(3) 



The retarded and advanced quantities, F = g, £<, II, are defined as 




(4) 




(5) 
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v: b R/A (i,i , )=v ab (i-i')6(t 1 -t[)+ 

+ ^Kc(i-i') J d2nf/ A (i,2)y^(2,i')- (6) 

cd 

(V a b denotes the Coulomb potential between charged particles of species "a" 
and "b"). Eq. (Q) has to be supplemented by its adjoint. Also, we will not 
consider correlated initial states here, for the corresponding generalization of 
the Kadanoff-Baym equations, see ||. 

To close this system of equations, we have to specify the carrier and longi- 
tudinal photon selfenergies £ and II. The simplest ansatz is the random phase 
approximation 

E|(l,lO=i^(l,l')^(l,lO, 

U*{l,l') = -ihg*{l,l')gf{l,l'). (7) 

Despite its simplicity, the system (||) - (Q) || describes remarkably complex 
physical processes: the evolution in space and time of charged carriers inter- 
acting via the full dynamic Coulomb potential which in turn evolves selfcon- 
sistently (screening build up) and may have nonequilibrium modes (including 
instabilities and nonlinear phenomena). Furthermore, the dynamics is influ- 
enced by the transverse electromagnetic field A which contains external fields 
(e.g. a laser field) and induced contributions and obeys Maxwell's equations. 
In this classical treatment, A does not modify the equations for carriers and 
plasmons explicitly, but only indirectly, via the particle propagators g R / A . 

Thus, a direct solution of the system (||) - (0) would yield a tremendous 
amount of information. However, here we proceed differently - we derive equa- 
tions of motion for the Wigner distributions of the charged particles f a (k, t) 
which turn out to be simpler. To this end we consider the correlation func- 
tions g< on the time diagonal. To make the derivations independent on a 
particular gauge, it is advantageous to make the Green's functions explicitly 
gauge-invariant |§,|o|. For the case of a spatially homogeneous field E(t) it is 
convenient to choose 4> ext = and A(t) = — c j t _ oo dt'E,(t). Then, the gauge- 
invariant Green's functions are obtained by the transform 



<7o(k, oj;t) = J drdr exp 




9a(r,r;t), (8) 



where r = r\ — r[, t — t± — t\ and t = (ti + t[)/2. The resulting quantum 
kinetic equation is given by [|10| 

d_ 

df 



r / Q (k ,t) + e a E(i) • V k / a (k a ,<) = -2Re f dt{z>g< - E<g>\ = J a , (9) 
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where the full arguments of the functions on the r.h.s. are, 

T$gf = T%[k* + Kt(t,fyt,i] gf[k a + Ki(t,i);i,t], (10) 

with K^^/V A(t) ' A(0 
a K ' c J t > t-t' 

Finally, to obtain a closed equation for the Wigner functions, the two-time 
functions and gj* in Eq. (g) have to be expressed in terms of f a . The 
appropriate solution is the generalized Kadanoff-Baym ansatz |fil"| which, in 
the case of time-dependent fields, has the form ]lO|JH 

± 5 |(k; tut',) = 3 f (k; hA) /I [k - Kftf, t); t[] 

-/| [k-Kf(t,0;ii] g^^tuA), (11) 

where the upper (lower) sign refers to g > (g < ), and f < = f and f > = 1 — /. 

Equations (^) and ([!(]) are valid for any choice of the selfenergies. In 
particular, we can use the RPA expressions (^), which leads to a collision 
integral I a || which generalizes the result of Haug and Ell JlJ,[b| to the case 
of an external electric field. A particularly interesting phenomenon is that the 
dielectric and screening properties of the plasma are directly modified by the 
electromagnetic field: 

n«M;t,t>) = - l -e(t-t>)e^^ J (*-*') 

x {f a [k + Q a (t,t');t'} - f a [k + q + Q a (t,t');t'}} , (12) 

where we used for g R / A in Eq. (11) the propagators of a free particle with energy 
e in an electromagnetic field [14|. Further, Q a and R a are, respectively, the 
momentum gain and displacement of a free particle in a field E(t) during the 
time interval [t' , t] given by 

Q a (t,t') = -e„ / dt"E(t"), R a {t,t') = / dt [ diE(t). (13) 

Jt' m a Jt' Jt' 

Obviously, the plasmon spectrum is modified in two ways: first, by the field- 
dependent prefactor of the standard RPA-polarization and, second, by the 
field-dependent momentum arguments of the distribution functions. 

For the practical issue of plasma heating it is important to understand 
the scattering process of charged particles in the presence of the field E(t) in 
detail. Analytical and numerical progress can be made if the RPA selfener- 
gies (0) are approximated by their static limit, which in Eq. (0) leads to the 
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substitution of the dynamic potential by the statically screened Debye poten- 
tial, Vf i ^ A (q,ti,t' 1 ) — > V st {q,ti)5{t\ — t[). As a result, the collision integral 
becomes Q 

I a (K,t) = 2j2 J ^^|Kf(k a -k Q )| 2 ( 27 r^(k a+ k b -k a -k b ) 
x J dt cos |^ [(e ab - e ab )(t - t) - (k a - k a )R 6(t,t)] | 

x {/.A [i - /«] [i - A] - /.A [i - /.] [i - /i]} | f , (14) 

where we denoted e a b = e Q + e&, e Q = p a /2m a , fa = f{K ai t), f a = f(K a ,t), 
K a = k a + Q a and R ab = R a — R b . Although plasmon and screening dy- 
namics are no longer included, Eq. ([)]) with the collision term (|lj) still con- 
tains important physics reflecting the influence of the electromagnetic field: 
(I) field-induced change of the arguments of the distribution functions, i.e. 
time-dependent generalization of the intra-collisional field effect; (II) modifi- 
cation of the energy balance in two-particle scattering [argument of the cosine 
in Eq. (|l4|)] - a strong field may essentially modify the energy broadening in 
the electron-ion integral; (III) nonlinear (exponential) dependence of the colli- 
sion integral on the field strength - this leads to the generation of higher field 
harmonics in the scattering processes. Furthermore, it gives rise to scatter- 
ing processes which involve emission (absorption) of photons, i.e. (inverse) 
bremsstrahlung. Indeed, it is straightforward to show |l0) that, for an electric 
field E(t) = Eocosflt transport quantities computed from ([hi]) will contain con- 
tributions proportional to J^(z)S[e a i, — l ao + qw„(,(i) — nMl], (— oo < n < oo), 
where the amplitude of an n-photon process is given by the Bessel functions 
J„. The argument z of J n is determined by the field strength and frequency, 
z = q[v° - vg]/fifl, where v° = e a E /m a tt and w afc = [v° - v£] sinOt. 



3 Numerical results 



The kinetic equation with the collision integral (14) is a convenient starting 
point for numerical investigations. We underline that it fully contains quantum 
effects, therefore, in contrast to classical equations, no cutting procedures at 
large momenta are required. Before presenting results of direct numerical 
solutions, we outline the idea for an approximate treatment |9|,|l5|]. 

If in Eq. (^J) the collision term is much smaller than the field term on the 
l.h.s., one can use an ansatz proposed by Silin |l6|], f a = /„+/„, where /° obeys 
the collisionless equation, with the solution /°(k, t) = f a o [k + ^■A(t)], and 
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the correction follows from Eq. (^|) with I a [f a ] — > I a [/°] . Further simplifica- 
tions are possible if the distributions are Maxwellians. Then, most of the inte- 
grations in the collision term can be performed, and analytical results for field- 
dependent transport quantities can be derived such as the conductivity, Fig. 1, 
or the cycle-averaged energy gain per particle, ( j • E ) = y J T dt' j(t') ■ E(i'), 
for which one obtains Pq| 

oo 

&V^Z 2 e A n e n l JmZ „ 2 2 f k f n 2 m e Q 2 

<J' E ) = ^3/9 n > " / rffc ■ ,.2» ex P- 



{k B Tf/ 2 ^ y {k 2 + k 2 D ) 2 *\ 2k B Tk 2 







X l-5SSf ) -gig 1 - I o* ■>•■ >) ■ < 15 > 

with the Debye screening wave number fc|j = A^ 2 = Airn e e 2 /ksT. We un- 
derline that quantum effects occur in ( |l5| ) in two places: first, the exponential 
function in the second line automatically ensures the convergence of the integral 
for large k and, second, the Bose-Einstein statistics of photons is reflected by 
the factor with the sinh function. For a classical plasma, both factors approach 
1, and one recovers the result of Klimontovich ]H| . Fig. 1 shows results for the 
conductivity computed using the Silin ansatz with the statically screened and 
the dynamically screened e-i collision integral. Also, we show a fit formula of 
Silin which holds for classical plasmas, but fails at low temperatures. 

Finally, we present results of direct numerical integration of the quantum 
kinetic equation (|^) with the collision term (|l4|). This has the advantage 
that no assumption on weak collisions or shape of the distributions has to be 
made. The heating of the electrons for three different densities is shown in 
Fig. 2, and Fig. 3 shows the evolution of the electron distribution during the 
first laser cycle. One clearly sees the anisotropy of the distribution and its 
broadening due to absorption of laser energy. A detailed analysis shows that, 
in fact, multi photon absorption (inverse bremsstrahlung) occurs which will be 
discussed elsewhere. 
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Figure 1: Conductivity of a hydrogen plasma with density n = 5- 10 21 cm 3 vs. temperature 
in a strong high frequency laser field with A = 248. 5nm, U/lu p i = 1.9, vg/vth = 5. 
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Figure 2: Electron heating due to Coulomb collisions in a strong laser field for different 
densities: n e = 10 22 cra~ 3 (full line), n e = 10 23 cm -3 (dashes) and n e = 10 24 cm~ 3 (dots). 
Remaining parameters same as in Fig. 3. (The straight lines are a guide for the eye.) 
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Figure 3: Evolution of the electron distribution during the first laser period. E(t = 0) = 
10 8 V/cm, figures correspond to subsequent quarter periods. Parameters are for hydrogen, 
see inset. 
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